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Atom loss resonances in ultracold trapped atoms have been observed at scattering lengths near 
atom-dimer resonances, at which Efimov trimers cross the atom-dimer threshold, and near two- 
dimer resonances, at which universal tetramers cross the dimer-dimer threshold. We propose a new 
mechanism for these loss resonances in a Bose-Einstein condensate of atoms. As the scattering length 
is ramped to the large final value at which the atom loss rate is measured, the time-dependent scat- 
tering length generates a small condensate of shallow dimers coherently from the atom condensate. 
The coexisting atom and dimer condensates can be described by a low-energy effective field theory 
with universal coefficients that are determined by matching exact results from few-body physics. 
The classical field equations for the atom and dimer condensates predict narrow enhancements in 
the atom loss rate near atom-dimer resonances and near two-dimer resonances due to inelastic dimer 
collisions. 
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Nonrelativistic particles whose scattering lengths are 
large compared to the range of their interactions ex- 
hibit universal low-energy behavior [T]. The universal 
few-body phenomena can include a spectrum of loosely- 
bound molecules as well as reaction rates of the particles 
and molecules. In some cases, including identical bosons, 
the universal behavior is governed by discrete scale invari- 
ance. The universal molecules then include sequences of 
3-particle clusters (Efimov trimers) [H [3] , 4-particle clus- 
ters (universal tetramers) [IHi], and clusters with even 
more particles [7] ■ The universal reaction rates exhibit 
intricate resonance and interference features [1]. 

The technology for trapping atoms and cooling them 
to extremely low temperatures has made the universal 
low-energy region experimentally accessible. The use of 
Feshbach resonances to control the scattering length a ex- 
perimentally makes ultracold atoms an ideal laboratory 
for universal physics. One particularly dramatic probe of 
universality is the loss rate of atoms from a trapping po- 
tential. Resonance and interference effects in few-body 
reaction rates can produce local maxima and minima in 
the loss rate as a function of a. The most dramatic sig- 
nature for a universal A^-atom cluster with A^ > 3 is 
the resonant enhancement of the A^-atom inelastic col- 
lision rate at a negative a where the cluster crosses the 
A^-atom threshold and becomes unbound. We refer to 
such a scattering length as an N-atom resonance. The 
first observations of an Efimov trimer [8 and a universal 
tetramer [9] and the first evidence for a universal 5-atom 
cluster [lU] were all obtained using a thermal gas of ^'^■^Cs 
atoms by tuning a to 3-atom, 4-atom, and 5-atom res- 
onances, respectively. An Efimov trimer has also been 
observed as an enhancement in the loss rate in a mixture 
of ^'^^Cs atoms and dimers [11] at a positive scattering 
length a* where an Efimov trimer crosses the atom-dimer 
threshold and becomes unbound. We refer to a* as an 



atom-dimer resonance. Another dramatic loss feature at 
positive a is an interference minimum in the 3-atom re- 
combination rate into the shallow dimer, which was also 
first observed in a thermal gas of -'^^■^Cs atoms [S]. Many 
of these loss features have been subsequently observed in 
ultracold trapped atoms of other elements [I2] . 

There are a few loss features in ultracold atoms that 
have not yet been related to universal few-body reaction 
rates. They all appear in systems that were believed to 
consist of atoms only and no dimers. Narrow enhance- 
ments of the loss rate near atom-dimer resonances have 
been observed in both a Bose-Einstein condensate (BEC) 
and a thermal gas of '^^K atoms [13] and in both a BEC 
and a thermal gas of ''Li atoms [TJl [TH]. In a BEC of 
^Li atoms, narrow enhancements of the loss rate have 
also been observed at positive values of a near two-dimer 
resonances [14] . at which universal tetramers cross the 
dimer-dimer threshold and become unbound. No mech- 
anism has been proposed for a narrow loss feature at a 
two-dimer resonance in a system consisting of atoms only 
and no dimers. One proposed mechanism for a narrow 
loss feature near an atom-dimer resonance a* in such a 
system is the avalanche mechanism, in which the 3-body 
recombination rate into the shallow dimer is amplified by 
secondary elastic collisions of the outgoing dimer [T^ . It 
was recently shown that the avalanche mechanism can- 
not produce a narrow loss feature near a* |16| . This 
is a consequence of the universal energy dependence of 
the atom-dimer cross section. Instead of having a nar- 
row peak at a* , the elastic cross section for the energetic 
dimer from the recombination event has a broad maxi- 
mum near 4.3 a*. Detailed Monte Carlo simulations of 
the avalanche process demonstrate that it does not pro- 
duce any narrow loss features [T7] . 

In this Letter, we propose a new mechanism for nar- 
row loss features near atom-dimer resonances and near 
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two-dimer resonances in a Bose-Einstein condensate of 
atoms. The mechanism is motivated by the phenomenon 
of atom-molecule coherence^ which involves the coherent 
transfer of atom pairs between an atom condensate and 
a coexisting dimer condensate. A small condensate of 
shallow dimers can be produced coherently by the time- 
dependent scattering length as it is ramped to the large 
final value where the atom loss rate is measured. The 
loss features then arise from the resonant enhancement of 
inelastic collisions involving dimers from the dimer con- 
densate. 

The phenomenon of atom-molecule coherence was dis- 
covered by Donley et al. in 2002 using a BEC of *^Rb 
atoms |18j . In these experiments, a pulse in the mag- 
netic field brought the atoms very close to a Feshbach 
resonance. The atoms were allowed to evolve at a large 
constant scattering length for a variable holding time, 
and then a second pulse took the atoms close to the Fes- 
hbach resonance again. Subsequent measurements of the 
number of atoms revealed three distinct components: a 
"remnant" BEC, a "burst" of relatively energetic atoms, 
and "missing" atoms that were not detected. The num- 
bers of remnant, burst, and missing atoms all varied si- 
nusoidally with the holding time at the frequency asso- 
ciated with the dimer binding energy. That sinusoidal 
dependence can be explained by a coexisting condensate 
of shallow dimers that was created by the first pulse [TW - 
[5T]. 

The behavior of atoms and dimers with sufficiently 
small kinetic and potential energies can be described by 
a quantum field theory with independent fields for the 
atoms and dimers. Atom and dimer condensates are 
described by classical fields tp{r,t) and d{r,t) that are 
the expectation values of the quantum fields. The quan- 
tum field equations can be formulated as coupled integro- 
differential equations for ip, d, and an infinite hierarchy of 
correlation functions for quantum fiuctuations. A typical 
experiment on the atom loss rate in a Bose-Einstein con- 
densate begins with a stable BEC of atoms with a small 
positive scattering length. This system is described by 
a static atom condensate ilj{r), with d{r) = and zero 
correlation functions. A ramp in the magnetic field pro- 
duces a time-dependent scattering length with a large 
final value a. During the ramp, the system is described 
by time-dependent condensates ipif^t) and d{r,t) and 
nonzero correlation functions. At the end of the ramp, 
the dimer condensate ci(r, 0) will be nonzero, although it 
could be very small if the ramp is nearly adiabatic. It 
could presumably be calculated using the methods de- 
veloped in Refs. P^U^ET] to describe atom-molecule co- 
herence. We will not attempt to calculate c?(r,0), but 
simply take the initial fraction fu of the atoms that are 
bound into dimers in the dimer condensate to be an un- 
known initial condition. We will ask whether observable 
loss features can be produced by the dimer condensate 
for a plausibly small fraction fu. 



During the subsequent holding time, there are tran- 
sient effects in the atom and dimer condensates and 
in the correlation functions that will die away. We 
will assume that after they have died away, the system 
can be described by coexisting atom and dimer conden- 
sates only and that the correlation functions can be ne- 
glected. Atom-molecule coherence will produce oscilla- 
tions in tp{r,t) and d{r,t) at the frequency fi/ {2'Kma?) 
associated with the dimer binding energy. The time- 
averaged condensates also decrease with time due to 
loss processes. The coexisting condensates can be de- 
scribed by a low-energy effective field theory for atoms 
and dimers whose kinetic and potential energies are small 
compared to the dimer binding energy. The interaction 
terms in the classical Hamiltonian density include 
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The coefficients of the interaction terms are univer- 
sal functions of a and the complex Efimov parameter 
K* exp(ir7H,/so), which can be interpreted as the bind- 
ing wavenumber of an Efimov trimer in the unitary limit 
a = ±oo [1]. The small positive parameter 77* takes into 
account tightly-bound diatomic molecules (deep dimers) , 
which provide decay channels for Efimov trimers. The 
coefficients of the interaction terms in Eq. ([T]) are con- 
strained by discrete scale invariance to be log-periodic 
functions of k*, with discrete scaling factor e^/'*", where 
So ~ 1.00624 is a universal constant. 

The coefficients in Eq. ([I]) can be determined by de- 
manding that few-body results in the fundamental theory 
be reproduced by the effective field theory. The coeffi- 
cient of d*d is determined by matching the rest energy 
of a shallow dimer: V2 = —Ti^ /ma?'. The coefficient of 
the terms that are 4"^ order in the fields can be de- 
termined by matching elastic scattering amplitudes at 
threshold: A2 = Stta, — STra^D, and fi = Airao, 
where a ad and au are the atom-dimer and dimer scat- 
tering lengths. The atom-dimer scattering length aAD 
is a multiplied by a simple log-periodic function of aw* 
[T]. In the limit ry* = 0, that function is real and it 
diverges at the atom-dimer resonance a* = 0.07076/kh,. 
The dimer scattering length ajj is a multiplied by a log- 
periodic function of ok* that is complex even if 77* = 0. 
For Ty* = 0, Deltuva has calculated a d with several digits 
of accuracy [22] . Its imaginary part has narrow resonant 
peaks at the two-dimer resonances Oi* « 2.196 0* and 
a2* ~ 6.785 a*. These two-dimer resonances were first 
calculated in Ref. [S]. For 77* > 0, a_D can be obtained 
by constructing an analytic fit to Deltuva's results for 
a£)/a as a function of ok* and then carrying out the an- 
alytic continuation k* — )■ k* exp(i77*/so). The coefficient 
of {ip*Tp)^ in Eq. ([ij could be determined by matching 
the elastic 3-atom scattering amplitude in the low-energy 
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limit. By the optical theorem, its imaginary part is pro- 
portional to the 3-atom recombination rate coefficient, 
ImAa = — (3m/?i)a, which can be separated into contri- 
butions from recombination into the shallow dimer and 
into deep dimers: a = ashaiiow + ctdoop- They both have 
the form ha'^/m multiplied by log-periodic functions of a 
[T]. In the limit 77* = 0, adocp = and ashaiiow has an 
interference zero at a+ — 0.31649/k*. 

The time dependence of ilj{r,t) and d{r,t) is deter- 
mined by the classical field equations associated with the 
interaction Hamiltonian density 'Hint in Eq. ([T]). The 
corresponding number densities are ua = and njj — 
d*d. The rates at which the numbers Na = Jd'^ruA 
and Njj = J d^rriQ change is determined by the anti- 
hermitian part of Hint: 
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Since unitarity requires the imaginary parts of uad and 
to be negative, these equations imply that Na and 
No both decrease monotonically with time. This ex- 
cludes atom-molecule coherence, which involves coherent 
oscillations between Na and Njj with angular frequency 
Ed jfi. The appropriate interpretation is that the high- 
frequency variations in the condensates associated with 
atom-molecule coherence are not resolved within our low- 
energy effective field theory. It can at best describe num- 
ber densities ?T.^(r, and 7i£)(r, t) that are averaged over 
many periods of the atom-dimer oscillations. 

To predict the loss rate of atoms during the holding 
time, we need initial conditions -0(r, 0) and d(r, 0) that 
are robust approximate solutions of the classical field 
equations associated with the hermitian part of Hint in 
Eq. ([!]). At the start of the holding time, the scattering 
length has its final value a and there is a specified initial 
total number of atoms A'o, an unknown fraction /^i of 
which are bound atoms forming the dimer condensate. 
We consider atoms trapped in a cylindrically symmet- 
ric harmonic potential: V{r) — ^mu!l[z'^ + + y^)]. 
If fo is sufficiently small, the effect of the dimer con- 
densate on the atoms is negligible. If the ramp to the 
final scattering length is slow enough, the atom conden- 
sate remains adiabatically in its ground state. We as- 
sume Nau is large enough that the kinetic energy of the 
atoms is small compared to their potential and interac- 
tion energies. The atom condensate then has the familiar 
Thomas-Fermi density profile: 

nA{r) ^ ^ ^2 niax{0,/i^ - y(r)}. (3) 
4:Trn a 

The chemical potential is determined by the num- 
ber Na = (1 — fD)No of unbound atoms: ^a = 
imw^a^flSC^A^Afl/oz]^/^, where = {fi/mujzY/'^ . 
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FIG. 1: (Color online) Three-atom recombination rate coef- 
ficient 1/3 as a function of the scattering length a. The data 
points were measured using a BEC of ^Li atoms [23] . The ar- 
rows point to the narrow loss features identified in Ref. [23j . 
The lower and upper parallel curves are the universal result 
for L3 (with a+ = 1402 ao and -q, = 0.038) and 9.2 Ls- The 
three vertical lines are the universal predictions for a,, ai», 
and 02* using a+ as input. The additional atom-dimer and 
dimer-dimer contributions in L"^ are included in the thinner 
solid (green) and dashed (red) lines, respectively. In the re- 
gions near a,, ai«, and 02*, the dimer fractions are fo = 10~^, 
3 X 10"^ and 6 x 10"^ 



We next consider the initial condition for d{r, 0). Since 
the initial number of dimers /£)iVo/2 is small compared 
to A''o, their self-interaction energy is negligible but the 
mean-field energy from the atom condensate can be im- 
portant. Assuming the dimer condensate produced by 
the ramp of the scattering length is in its ground state, 
it satisfies the Schroedinger equation 
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The eigenvalue /iu is the chemical potential of the dimer. 
The normalization of the dimer condensate is determined 
by the condition / d^r d*d = fDNo/2. The total potential 
energy of the dimer is the sum of the trapping potential 
2V{r) and the mean-field energy. Its minimum is at the 
origin if KeaAD < fa and near the edge of the atom 
condensate if KeuAD > |a. In these two cases, the dimer 
condensate is an ellipsoid centered at the origin and an 
ellipsoidal shell near the edge of the atom condensate, 
respectively. The boundaries in a for the ellipsoidal shell 
region are at a* and 2.86 a*. 

The data on the loss rate of ^Li atoms in Ref. [T3] 
has recently been reanalyzed using a more accurate de- 
termination of the parameters of the Feshbach resonance 
[23] . Under the assumptions that the system is a pure 
BEC of atoms and that the loss rate comes from 3-atom 
recombination only, the rate coefficient L3 is defined by 
dN/dt = -(L3/6) / d^rn\. The data from Ref. I23j that 
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was measured in a BEC of ''Li atoms is shown in Fig. [T] 
The initial number of atoms was A'o = 4 x 10^. Their fit 
to the universal prediction L3 = 3a using an adjustable 
normalization factor is shown in Fig. [1] It deter mines 
the Efimov parameters a+ — 1402 qq and r?* — 0.038 
[13] and requires a normalization factor of approximately 
9.2. The corresponding universal predictions for the reso- 
nances are a* = 313 oq, ai» = 688 ao, and 02* = 2127 ao- 
Near each of these three resonances, there is a narrow loss 
feature where the data are significantly higher than the 
fitted curve. The positions of the loss features reported 
in Ref. are 426 ao, 919 ao, and 1902 ao. As shown 
in Fig. [1] they are not as widely spaced as the predicted 
resonances. 

We now consider the effects of a small coexisting dimer 
condensate. Our initial conditions are determined by 
iVo = 4 X 10^ and an assumed dimer fraction fjy, which 
we expect to be small compared to 1 and to depend 
strongly on a. We take ua to be the Thomas-Fermi 
profile in Eq. ([s]). For no = d*d, we use a varia- 
tional approximation for d that reduces the 3-dimensional 
Schroedinger equation in Eq. Q to a 1-dimensional equa- 
tion for a single radial variable. The subsequent time 
dependence of the total number of atoms N{t) can be 
obtained by solving Eqs. ([2]). A quantitative compar- 
ison with the data would require comparing with the 
time dependence observed in the experiment. Instead 
we compare the data for L3 in Fig. [T] with an effective 
rate coefficient determined by the initial loss rate. 
It is defined by dN/dt = -{Lf/6) Jd^r 

'^AO' where the 

integral is evaluated under the assumption that fo = 0: 
Jd^rn\^ = iNQ/168Tr^ata^)[15CNoa/a,]^/^. In Fig. |l| 
the universal prediction for L3 from 3-atom recombina- 
tion falls an order of magnitude below the data. This 
allows room for additional contributions from the atom- 
dimer and dimer-dimer terms in Eqs. ([2|. We would like 
to determine whether narrow loss features can stand out 
above the smooth contributions for plausible values of 

fD- 

In Figjl] the terms in L^^ are illustrated by two thin 
curves, which correspond to 3-atom recombination plus 
atom-dimer losses and to 3-atom recombination plus 
dimer-dimer losses. We show these curves near a«, ai*, 
and a2* using different values of fo chosen to make the 
narrow loss feature visible: fa — 10"^, 3 x 10"'^, and 
6 X 10~^, respectively. The curve near a* that includes 
the atom-dimer terms in Eqs. ([2]) has a narrow peak at 
the atom-dimer resonance. The peak is not symmetric, 
because the dimer condensate changes from an ellipsoid 
centered at the origin for a < a* to an ellipsoidal shell 
near the edge of the atom condensate for a > a* . The 
curves near ai* and a2* that include the dimer-dimer 
term in Eqs. ([2| have narrow peaks at the two-dimer 
resonances. The behavior near ai* and a2* is different, 
because near ai* the dimer condensate is an ellipsoidal 
shell while near a2* it is an ellipsoid. If fu is large enough 



near ai* and a2*, the peaks can stand out above the 3- 
atom recombination and atom-dimer contributions. If fo 
is too large near a+ , the atom-dimer contribution can fill 
in the interference minimum from 3-atom recombination. 
For r]^ = 0.038, there is no longer a local minimum near 
a+ if /d > 2.2 X 10^^. The dimer fractions illustrated 
in Fig. [1] are sufficient to make the atom-dimer and two- 
dimer resonances stand out in the initial loss rate. Larger 
fractions would be required to make a significant differ- 
ence in the integrated loss rates. Nevertheless our results 
demonstrate that narrow loss features can be produced at 
the atom-dimer and dimer-dimer resonances with plausi- 
bly small values of the dimer fraction. 

Our dimer condensate mechanism provides a plausible 
explanation for the narrow loss feature at an atom-dimer 
resonance that was observed in a BEC of '^^K atoms [T3] 
and for the narrow loss features near an atom-dimer res- 
onance and near two two-dimer resonances that were ob- 
served in a BEC of ^Li atoms [14]. It cannot explain 
the narrow loss features near atom-dimer resonances that 
have been observed in thermal gases of '^^K atoms [15] 
and ''Li atoms [T31 [H]. These loss features were ob- 
served at relatively small scattering lengths, so they could 
be associated with nonuniversal effects. In a BEC, the 
dimer fraction fjj would depend sensitively on the de- 
tailed form of the ramp that brings the scattering length 
to its final value a. This sensitivity could be exploited to 
test our mechanism. The fraction fu could be amplified 
by pulsing the magnetic field very close to the Feshbach 
resonance before measuring the loss rate, as in the experi- 
ments on atom-molecule coherence [18] . The loss rates at 
atom-dimer and dimer-dimer resonances would increase 
linearly and quadratically with the amplification factor, 
respectively. 

In summary, we have proposed a new mechanism for 
narrow atom loss features at atom-dimer resonances and 
at two-dimer resonances in a BEC of ultracold atoms. 
The positions of these features are determined by univer- 
sal few-body physics and thus provide additional tests of 
universality. There could be similar loss features where 
universal 5-atom clusters cross the atom-dimer-dimer 
threshold. The strengths of all these loss features are 
determined by the many-body physics of Bose-Einstein 
condensates and open a new window into the remarkable 
phenomenon of atom-molecule coherence. 

We thank R. Hulet for valuable discussions. This re- 
search was supported in part by a joint grant from the 
Army Research Office and the Air Force Office of Scien- 
tific Research. 
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